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Abstract 

A model which combines the perturbative behavior of QCD with low energy 
phenomenology in a unified framework is developed. This is achieved by ap- 
plying a similarity transformation to the QCD Hamiltonian which removes 
interactions between the ultraviolet cutoff and an arbitrary lower scale. Iter- 
ation then yields a renormalization group improved effective Hamiltonian at 
the hadronic energy scale. The procedure preserves the standard ultraviolet 
behavior of QCD. Furthermore, the Hamiltonian evolves smoothly to a phe- 
nomenological low energy behavior below the hadronic scale. This method has 
the benefit of allowing radiative corrections to be directly incorporated into 
nonperturbative many-body techniques. It is applied to Coulomb gauge QCD 
supplemented with a low energy linear confinement interaction. A nontrivial 
vacuum is included in the analysis via a Bogoliubov-Valatin transformation. 
Finally, the formalism is applied to the vacuum gap equation, the quark con- 
densate, and the dynamical quark mass. 
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I. INTRODUCTION 



There is currently an unfortunate dichotomy in analytical approaches to QCD. One either 
calculates using Lagrangian methods and, with the exception of Monte Carlo simulations, is 
limited to perturbative calculations, or one employs low energy models of QCD and eschews 
the known high energy behavior of the theory. There are only a few approaches where, 
for example, the calculation of hadronic form factors may be connected to the expected 
quark counting behavior at high energies. These include QCD sum rules NRQCD [Q, 
and to some extent light cone quantization This paper develops a formalism by which 
the known high energy behavior of QCD (pQCD) may be rigorously joined to standard 
phenomenological models of hadrons. This is achieved by using the renormalization group 
to run the scale of the QCD Hamiltonian to the hadronic regime where it smoothly joins a 
fixed phenomenological behavior. 

This approach imposes important constraints on the low energy portion of the model. 
If one wishes to recover pQCD, one must work with current quarks. Thus the low energy 
theory must also deal with current quarks. However phenomenology indicates that con- 
stituent quarks are the relevant degrees of freedom at low energy. Thus, one must allow for 
spontaneous chiral symmetry breaking in the model. This, in turn, means that a nontrivial 
vacuum must be incorporated into the theory. Finally, we see that one must use many-body 
techniques when calculating observables. The Hamiltonian-based renormalization group we 
employ is particularly useful in this regard because it allows for the use of nonperturbative 
many-body techniques. 

We call this approach the "Dynamical Quark Model" (DQM). The fact that the low 
energy behavior of the DQM evolves directly from pQCD is very helpful when analyzing the 
dynamical structure of the low energy portion of the model. For example, we are able to 
rigorously establish that the non-Abelian Coulomb interaction is responsible for the quark 
structure of confinement in the heavy quark limit. We are also able to resolve ambiguities 
in the separation of short and long distance phenomena when evaluating current matrix 
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elements. 

This brings us to a final important reason for employing the renormalization group im- 
proved (RGI) Hamiltonian. As stated, we are obliged to use many-body methods while 
working with the model. The major tools at our disposal are the Tamm-Dancoff approxima- 
tion (TDA) and the random phase approximation (RPA). These build, for example, meson 
bound states of the type qq or qq + qq{qqqq) + ■ ■ ■ respectively. If one wants to allow for 
mixing to, for example, qqg or qqqq states then these must be explicitly included in the 
calculation. Such coupled channel problems are rather difficult to solve. Thus it is im- 
perative to include as much of the physics of Fock sector mixing as possible directly into 
the Hamiltonian. For example, the phenomenologically important hyperfine splittings can 
be incorporated by iterating the qqg vertex, thereby adding an effective transverse gluon 
exchange operator to the Hamiltonian. Using the RGI Hamiltonian automatically achieves 
this goal. 

The remainder of this paper is organized as follows. A model based on the Coulomb gauge 
QCD Hamiltonian is introduced in Section II and its high energy behavior is analyzed. To 
do this we impose a chiral invariant regulator and then calculate the RGI Hamiltonian. 
The counter-term structure and possible perturbative and nonperturbative renormalization 
schemes are then discussed. In Section III we show how effective operators are constructed 
and examine the ultraviolet behavior of the quark creation operator, the axial current, 
and the scalar quark density. Section IV presents nonperturbative, one-loop, many-body 
calculations of the vacuum gap equation, the quark condensate, and the dynamical quark 
mass. We summarize and conclude in Section V. 

II. THE EFFECTIVE HAMILTONIAN 

Our starting point will be the Coulomb gauge Hamiltonian of QCD. There are many 
reasons for choosing to work in the Coulomb gauge. For example, as stated in the Introduc- 
tion, it is expedient to have explicit qq interactions in the Hamiltonian to make a connection 
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with quark model phenomenology. Further advantages of the Coulomb gauge are that all 
unphysical degrees of freedom have been eliminated so that subsidiary conditions on states 
are not necessary. This also implies that the norm is positive definite. Furthermore, spurious 
retardation effects are minimized. 

We note that it is impossible to address the issue of gauge invariance in this model 
because it is defined in a specific gauge. However, QCD is gauge invariant so it is perfectly 
acceptable to work in a fixed gauge. The real issue is not how the model changes under a 
gauge transformation, but how well the phenomenological low energy interaction Vconf we 
will use mimics the actual low energy behavior of QCD. The only way to address this is to 
solve QCD (on the lattice, say) or to compare the predictions of the model to observables. 
This will be the subject of future papers. For now, we note that lattice gauge theory indicates 
that Vconf should be a linear potential when the quarks are static, our predicted glueball 
spectrum agrees well with lattice calculations [Q, and the model is guaranteed to reproduce 
the success of the naive quark model for heavy quarks. 

The DQM Hamiltonian is obtained after a rather long series of operations. The first step 
is to regulate the bare Hamiltonian in a manner which is consistent with the renormalization 
scheme we shall employ and which does not interfere with the nonperturbative calculations 
we will eventually perform. This introduces a scale into the Hamiltonian which we shall 
call Aq. Renormalization is carried out by performing a similarity transformation on the 
Hamiltonian. This transformation is designed to eliminate couplings between the scale Aq 
and an arbitrary, but nearby, lower scale Ai. The similarity transformation preserves the 
spectrum of the Hamiltonian (to the order we work to) while restricting the matrix elements 
of the Hamiltonian to a smaller, band-diagonal form. Since elimination of these couplings 
reveals UV divergences in the bare Hamiltonian the third step is to absorb these by fixing 
appropriate counter-terms. After this we repeat the last two steps above, thereby generating 
a renormalization group improved effective Hamiltonian. The final step of our program is to 
include a phenomenological term which is meant to model the behavior of the higher order 
terms that have been ignored in the perturbative analysis. 
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A. Model Definition 



As mentioned above, our starting point is the Coulomb gauge Hamiltonian of QCD 



H = Ho + V, 



(1) 



where Hn is the free Hamiltonian, 



Ho= dx^lj{x.y [-ia ■ V + /Sm] ijj{x.) +Tr dx 



n2(x) + B^(x) 



(2) 



The degrees of freedom are the transverse gluon field A = A°T°, its conjugate momentum 
n, and the quark field in Coulomb gauge. We have represented the Abelian portion of the 
non-Abelian magnetic field by B^. The interactions are given by 



V = Tr 



+ Tr d:K 



B^(x)-B^(x) 



ju{^)j-'u{^) - n2(x) 

- g J d^ip\x.)a ■ A(x)^/'(x) (3) 

where = Det[V-D], D = V — (^A, is the covariant derivative in the adjoint representation 
and B = B, = VjAk - VkA^ + g[A^, A^]. 

The density, p, in the non-Abelian Coulomb interaction is the full QCD color charge and 
thus has both quark and gluonic components. 



p»(x) = ^t(x)TXx) + /"^'^A''(x) ■ n'^(x) 



(4) 



B. The Similarity Renormalization Scheme 

An effective Hamiltonian is defined to be an operator which, when acting within a sub- 
space of the original Hilbert space, yields the same observables as the original Hamiltonian. 
Usually the effective Hamiltonian is constructed by dividing the Hilbert space into two por- 
tions which are defined by projection operators, P, and Q = 1 — P. In general the effective 
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Hamiltonian may be written as a sum of the projected Hamiltonian, PHP and a set of ef- 
fective interactions, V^ff arising from elimination of couplings between the P and Q spaces. 
Typical methods for constructing the effective Hamiltonian produce energy denominators in 
Veff which may vanish. This introduces severe complications when attempting to nonpertur- 
batively diagonalize the effective Hamiltonian. An elegant way for avoiding these problems 
is given by the similarity transformation scheme of Glazek and Wilson 

The scheme involves ordering the states according to their free energy levels. The in- 
teraction part of the Hamiltonian is then regulated by ignoring matrix elements with free 
energy differences which are greater than Aq. One then proceeds by constructing a similar- 
ity transformation which removes the couplings between states whose energy differences lie 
between Aq and Ai. This procedure generates an effective potential which incorporates the 
eliminated physics. 

We start by constructing the eigenstates of the free Hamiltonian, Hq, by expanding the 
quark and gluon field operators in normal modes, 

i'i^) = j:lj^,{u(k,TMk,T)+v{-k,T)d*{-k.T))e->'- (5) 

The quark operators are given in the helicity basis and all discrete quantum numbers (he- 
licity, color and flavor, for the quarks and color for the gluons) are collectively denoted by 
r and a respectively. Note the use of non-relativistic normalization so that, u^u = v'^v = 1, 
and 

{6(ki, n), b\k2, T2)} = {rf(-ki, n), d\-k2, T2)} = {27lfS{ki - k2)Sr„r, (7) 

and 

[a(ki,ai),at(k2,a2)] = {2n)^6{k, - k2) (l - kiki) 5a,a,. (8) 
The free Hamiltonian is then given by, 
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^o = E/ 7S3^(k) [&nk,r)6(k,r)+cit(k,rKk,r)] J ^.;(k)at(k, a)a(k, a 



dk 

(9) 



where E(k) = + k^, ti;(k) = |k| and we have dropped the constant, zero point energy. 

In the basis of eigenstates of H^, Ho\n) = En\n), the cut-off method of Glazek and 
Wilson may be written as 

{m\H\n) ^ (m|if^°|n) = EJ^m + fmniAo){m\V\n). (10) 

Here /mn (Aq) is a function which is unity for \Emn\ << Aq and zero for \Emn\ >> Aq (we 
have defined Emn = Em — En). For our purposes we find it convenient to use a sharp cut-off, 

/,n„(Ao)=^(Ao-|^„„|). (11) 

In this paper we will study only the quark sector of the Hamiltonian. Thus all gluon 
interactions will be ignored except those which contribute to effective quark operators. A 
discussion of the gluonic sector has been given in Ref. 0], where it has been successfully 
applied to the glueball spectrum. The effort in deriving effective operators in the gluon 
sector is currently in progress 0]. 

The operators of the bare regulated Hamiltonian contributing to order in the quark 
sector are, 

= gV^^" + g%^° (12) 

where 

_ ^ r dk, dk2 dl 



TiT2a 



'J (2vr)3 (27r)3 (2vr)3 ^ ^ 
[u\au2] e{Ao - \E{k2) + u;{l) - E{k,)\) b\k,, n)b{k2, r^Ml, a) + 



(13) 



and 
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X 



^(Ao - |^(k3) + E{k^) - E{k,) - E{k2)\) [-uU3][4u4 



6t(ki, ri)6t(k2, r2)6(k3, r3)6(k4, n) + ■ ■ 



(14) 



In Vi and V2 the ellipses denote an additional seven and fifteen terms respectively and Vc is 
the momentum space Coulomb potential. The additional terms correspond to interactions 
involving all possible permutations of quark and antiquark operators. The subscripts on 
the spinors and the Gell-Mann matrices refer to indices of the appropriate quark operators. 
Normal ordering the bare Hamiltonian with respect to the perturbative basis generates a 
new one body operator given by 

r dk r 

g'Vseif.c = E y 7^ [Sc(k; A?)[&t(k, T)b{k, r) + ^;t(-k, r)rf(-k, r)] 



- E{k))Gcik- Ao^) [6t(k, r)dt(-k, r) + d(-k, r)6(k, r)]] . (15) 



Both S(7 and Gc (which are given in the Appendix) are divergent and not regulated by Aq. 
In Eq. [^] they have have been regulated individually with a cutoff denoted by A^. Since all 
divergences of the bare Hamiltonian are to be absorbed by counter-terms one could subtract 
this self energy in its entirety and simply ignore it. However we shall show subsequently 
that covariance requires that a portion of KeZ/.c be maintained in the Hamiltonian. The 
Oi^g^) bare potential will therefore be given by, 

r^"(A^) = gV,^^^ + + /Kel/.clAo"^) (16) 

where we have made the Aq and A^ dependence exphcit. The Aq dependence of Vgeif.c arises 
only from the term proportional to Gq since Sc is diagonal in the perturbative basis and 
therefore, according to Eq. [|T0|, is Aq independent. So long as the ratio A^/Aq remains 
finite as both scales approach infinity, the existence of a second scale will be irrelevant (to 
the order to which we work). From now on we will set Af = Aq. 

We are now ready to calculate the equivalent of the P-space Hamiltonian, PHP = H^^ 
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and of the effective Hamiltonian, H^jf- The matrix elements of the former are defined using 

Eq. irg, 

(m|i/^» = EJnm + ^n(Ai)(m|y|n), (17) 

where V depends on Aq (and in general) as discussed above. Since within the P-space the 
spectrum of the bare Hamiltonian and that of the effective Hamiltonian are to be identical, 
we may write 

H^}f = SH^'^S-^ (18) 

where is a unitary matrix such that matrix elements, {m\H^jj-\n) are nonzero only within 
the P-space, i.e. for states that satisfy 

\Enm\ < Ai. (19) 

To construct S we proceed by expanding it in powers of the strong coupling 

S = e'^, (20) 

and 

R = gRi + g^R2 + 0{g'). (21) 

Then 

+ {V2 + Vseif.c + i[R2, Ho] + z[Pi, V,] - i[Pi, [Pi, Po]]) + 0(^7=^). (22) 

Here, the dependence of the V^s on Aq is implicit. To order g, ig[Ri,HQ\ is chosen to 
eliminate matrix elements of gVi between the scales Ai and Aq. Thus, to this order, one 
need simply replace the two terms in the first brackets with V^^ , 

{m\V,^'\n) = e{A, - \E^n\){m\V,\n). (23) 
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Similarly, to order (7^, i(7^[i?2, -f^o] eliminates all matrix elements of the expression in the 
second brackets between Aq and Ai. A new interaction in H^jj which distinguishes it from 
H^^ arises from the double commutator term in Eq. ||2^ when it is evaluated for matrix 
elements satisfying Eq. H^jj is thus given by 

Htif = Ho + gV,^^ + g^Vt + g^V,%,c + /^t ^ + 0{g^). (24) 

where 

{m\Vt + Ke//.cl^) = - \E^n\){m\V2 + Kei/.c|n). (25) 

The last term in Eq. [24] has been generated by the similarity transformation and is meant 
to reproduce the physics lost when the cutoff is reduced from Aq to Ai. For sharp cutoffs 
(see Eq. |]lT|) one has 

{m\V^'\n) = e{Ai - |E„„|)^(m|ri|g)(g|ri|n)eT(^,n,^,rr.;Ai). (26) 

Note that the sum over intermediate states is cut-off from above due to an implicit depen- 
dence of the matrix elements of Vi on Aq. The last factor is a function which arises in the 
similarity transformation and is similar to an energy denominator in standard perturbation 
theory. It is given by 

a,; A) = (^1 _ 1 A,| - A)) + (a ^ 6) (27) 

Notice, that energy denominators never cause divergences since they are cut off from below 
by Ai. 

We are now in a position to study the ultraviolet behavior of the theory. As Aq approaches 
infinity V^^ is finite since it is cutoff at the lower scale Ai. Alternatively, the self energy 
Vseif.c diverges as A^ = Aq ^ 00. Similarly, the effective interaction V^^ contains self 
energy pieces which diverge and interaction pieces which are ultraviolet finite. Thus it is 
convenient to split it into two portions 

^ = Kt/V.T + Vtr, (28) 
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where 

9X%.T = i:J ^ET(k;Ao,Ai)[6nk,r)6t(k,r) + dt(-k,r)d(-k,r)] 

+ - E{k))GT{k; Ao, Ai)[b^ik, T)d^{-k, r) + d{-k, r)6(k, r)] (29) 

and Vex.T has a similar structure to V2. The functions Et and Gt are given in the Appendix. 
As stated above, in the hmit Aq ^ 00, V^lj^ j, is divergent while the exchange interaction, 
Kt^T remains finite. 

The divergences in the self energy pieces must be absorbed by adding appropriate 
counter-terms to the bare Hamiltonian. Thus we define 

H = lim \h^° + g^Vj^°] . (30) 

Since only one-body operators are explicitly UV divergent we construct counter-terms of the 
form 

gXt-T.] ^S^*-(k;Ao)[6t(k,r)6(k,r)+dt(-k,r)d(-k,r)] 

+ - E{k))G'^-{k; Ao) [b\k, r)dt(-k, r) + d{-k, r)6(k, r)]. (31) 

Note that Vj^° contains a piece which is off-diagonal in the perturbative basis {G'^-). This 
will lead to a contribution from in R2 when -|- g'^V^° is used as an input to the 
similarity transformation (as it should be). This extra contribution to R2 removes couplings 
caused by G'^^- when Ai < 2E{k) < Aq and therefore only leaves a nonvanishing contribution 
from this counter-term for Ai > 2E{k). 

C. The Renormalization Group Improved Hamiltonian 

Our goal is to define an effective Hamiltonian which is operative at hadronic scales. Thus 
one may be tempted to set Ai ~ 1 fm~^. However, a single-step similarity transformation 
which evolves the Hamiltonian from the UV cut-off Aq down to the hadronic scale would in- 
troduce effective interactions with large matrix elements for which the perturbative approach 
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is hard to justify. This would also result in a large sensitivity of the effective Hamiltonian 
to parameters which are set at the UV scale. Instead we may use the renormalization trans- 
formation iteratively over many small steps and at each step adjust the Hamiltonian only 
by a small amount. In the language of standard Feynman perturbation theory, the two 
approaches correspond to keeping only the leading-log and summing the leading-log to all 
orders, respectively. Thus we choose to apply the similarity transformation 



logfl + e) 



(32) 



times starting at the beginning of the renormalization group trajectory with the bare Hamil- 
tonian, = H^^jr and finishing at the end of the trajectory with H^^'j. Each time the 
transformation is applied the cutoff is reduced from Aj to Aj+i by a finite amount 

A. 



A, 



1 + e. 



(33) 



The limit Aq ^ oo corresponds to A^e diverging logarithmically. 

Expanding the integrands in the expressions for KeZ/.c and Ke«/.T (Eqs. [15] and [29]) 
in powers of e yields (see Appendix), 



Gc(k;Ao 



Er(k; Ao,Ai) 
Sc(k; Ao) 
Gr(k;Ao,Ai) 
Gc(k;Ao,Ai) 



2m2 



8k2 



E{k) 3E{k) 

4m2 8k2 

+ 



1 + 



1 



EiX) 3E(k) 



l + 0(^))+Sc;(k;Ai) 



^(47r)2 

g'^Cp 14|k|m ^ / Ai \ A ^ ^ / 1 

o^C^ 4|k|m^ /Ai ^„ ,\ / ^fl 
+G'c;(k;Ai,Ai). 



(34) 



The presence of the last terms in the second and fourth equations above causes a complica- 
tion. Since A^r is fixed at the end of the trajectory these terms are of 0{Ne) and dominate 
transverse self energies generated by the similarity transformation, which are of 0(e) (the 
first and third equations of Eq. |p3). In order for these two types of self energies to be 
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comparable, most of the Coulomb self energy induced by normal ordering must be sub- 
tracted by a counter-term in the original Hamiltonian. An unusual feature of this is that the 
subtraction must depend on the renormalization group trajectory (as opposed to standard 
perturbation theory). In other words it is not enough to fix the Coulomb counter-term in 
if'^o, it must also depend on the trajectory {i.e., it must be e-dependent). This is to be 
expected: covariance relates both Coulomb and transverse gluon exchange and since the 
similarity transformation affects only the latter, Coulomb self-energies must be adjusted in 
an appropriate fashion along the entire trajectory. 

For a Hamiltonian evaluated at Aj along a trajectory defined by Eq. |^3|] the counter-term 
is given by 



The first portion of each equation is meant to indicate that the counter-terms have become 
trajectory dependent. Possible finite pieces have not been considered in these expressions 
{i.e., scheme dependence) since they become irrelevant upon use of the renormalization 
group equations. 

At this stage the perturbative portion of the Hamiltonian is fully specified and we are 
ready to apply the renormalization group procedure. After a single step (Aq ^ Ai) a new 
one-body operator is generated. 




(35) 



and 




(36) 




(37) 



We have defined new quantities 



E(k;Ao,e 



) = Sr(k;Ao,Ai) + Sc(k,Ao) + 



S^*-(k,Ao,e) = e 



(47r)2 E{k) 



(H-O(^)) (38) 
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and 



^(k; Ao, e) = Grik; Aq, Ai) + GciK Aq) + G^* (k, Aq, e 



— e 



g'^Cp 6|k|m 
(47r)2 E(k) 



1 + 



(39) 



The last equality in each equation follows from Eqs. [34], [35], and [36]. In order to iterate 
the similarity transformation further the one body operator (Eq. [37]) must be diagonalized. 
This is because the eigenstates of this operator must be used as a new basis in which matrix 
elements of the interaction are regulated. We therefore proceed by rotating the quarks back 



to the massless basis (6, d ^ b,d). Eq. then becomes 



dk 

—mi 



\k\ + 0{- 
dk ^,Ai 



[b\k, T)b(k, r) + S{-k, T)d{-k, r)] 



(27r) 



E(k))[6t(k, T)d^ i-k, r) + d{-k, r)6(k, ■ 



(40) 



where 



mi = m 



6Cj 



(41) 



Note that it is only necessary to rotate the basis below the scale Ai so that the ^-functions 
may be ignored in this process. 

To order this looks just like the free Hamiltonian with m replaced by mi and quantized 
in a massless basis. We now rotate the one body operator back to the massive quark basis 
using mi as the bare quark mass and obtain a free, diagonal Hamiltonian, H^jj{mi) which 
may be used as input to the next similarity transformation (Ai A2). Note that since 
the similarity transformation is only affected to order g"^ by these rotations, the net effect 
of this procedure is simply to replace mo = m(Ao) by mi = m(Ai). The remainder of the 
renormalization group procedure is simple: we apply the above transformations times and 
obtain an effective Hamiltonian with the one-body operator given by the free Hamiltonian 
with the replacement mo — mAr = m(A7v): 
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mo — >■ rriN = m 1 + g 



2 6Cf 



(42) 



(47r) 



In the limit Ne ~ logAo/Ajv — >^ oo this becomes 



m(AAr) = m(Ao) (-r-^ 




(43) 



Finally, the remainder of the RGI Hamiltonian (the two-body terms) is given by the 
interactions defined at the end scale 



It is perhaps useful to summarize the calculation to this point. We have started with 
the exact Coulomb gauge Hamiltonian of QCD, expanded it to order g'^, and regulated it 
in a manner which lends itself to the similarity transformation renormalization procedure. 
Performing the transformation generated Coulomb and transverse gluon exchange self energy 
effective interactions which diverge in the ultraviolet limit. The divergences were absorbed 
into counter-terms and the whole process was iterated to construct the RGI Hamiltonian. 
This Hamiltonian is defined at a hadronic scale A^v, is ultraviolet and infrared finite, and 
has nonzero matrix elements only for states which differ in energy by up to A^r (higher 
energy physics is incorporated in the effective interaction, Vex.r)- Notice that all self energy 
interactions have been absorbed by the renormalization procedure and no longer explicitly 
appear in the Hamiltonian. As expected for the Coulomb gauge, there is no need for wave 
function renormalization and the coupling constant does not run until 0{g^). 



The RGI Hamiltonian which we have constructed can be rehably employed to calculate 
observables down to a scale of A ~ lfm~^ where the coupling constant becomes large (note 
that we shall refer to the end point of the renormalization group trajectory as A from 
now on). At this point the Hamiltonian must both be solved nonperturbatively and must 



gV,^- + g'V,^- + g'V,% 



(44) 



D. The Phenomenological Interaction 
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be supplemented with a phenomenological potential which incorporates as much of the 
neglected higher order physics as possible. The latter forms the topic of this section. 

Choosing an effective low energy interaction is not trivial because the model is relativistic. 
This means that a Dirac structure must be imposed on the interaction - with commensurate 
implications on phenomenology. It is commonly held that the effective potential should be 
linear with a scalar current interaction. However, we have shown that this simple expectation 
is misleading One may establish the Dirac nature of the low energy effective potential in 
the heavy quark limit. In this case it is appropriate to make a nonrelativistic reduction and 
use the quenched approximation. P-wave splittings in the J/ip and T spectra imply that 
long range spin-dependent potentials correspond to a nonrelativistic reduction of a scalar, 
phenomenological, confining interaction This view is also supported by lattice gauge 
calculations |^ and phenomenological, minimal-area-law model for the Wilson loop 



Starting with the Coulomb gauge Hamiltonian, it is possible to show that the effective 
scalar interaction arises dynamically through mixing with hybrids while the static part of 
the qq potential remains of a timelike- vector nature 

If one assumes that the heavy quark potential also applies for light quarks then the 
net result is to simply replace the order Coulomb potential with a Coulomb plus linear 
potential: 

- 9^Cf^ ^ -9^Cf^ + (45) 
47rr 47rr 

We shall fix the string tension at a = 0.18 GeV^, commensurate with Regge phenomenology, 
lattice gauge theory, and the constituent quark model. Note that the phenomenological two- 
body potential is meant to reproduce low energy physics, thus Eq. is used in the effective 
Hamiltonian, i.e., below the scale A. We also introduce a one-body quark operator which 
corresponds to the linear term in Eq. [45]. This eliminates color non-singlet states from the 
spectrum. Furthermore, it ensures that color-singlet bound states are infrared-finite. 

These observations will be discussed in greater detail below. The final form for the quark 
sector of the model Hamiltonian which we shall use in subsequent calculations is: 
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HdQM = //o(A) + + ^V^^ + ^V,ir + Vconf + 1^.e//.con/. (46) 

This has been defined at the scale A; hence the interaction terms are all cutoff at this scale. 
As discussed above, the A dependence of Hq comes from the running mass. Finally Vconf 
has only two-body matrix elements. 



III. EFFECTIVE OPERATORS 

If one wishes to make statements about matrix elements which are accurate to order it 
is not sufficient to know the states (z.e., the Hamiltonian) to this order - one must also evolve 
the relevant operator. This is illustrated below with explicit calculations of the perturbative 
evolution of the single quark operator, the axial charge operator, and the composite operator 
qq used to define the chiral quark condensate. 

Just as with the Hamiltonian, one may use the similarity renormalization scheme to 
calculate effective operators restricted to P-space. Following the steps discussed in Sec. II 
we first define an operator O^" by restricting its matrix elements in the basis of cigcnstatcs 
of //o, 

{m\0\n) (m|C^°|n) = ^(Aq - \Emn\){m\0\n) (47) 

The effective operator is then given by 

Oe// = SOS-' = O + ig[R,, O] + ([R^, O] - O]]) + 0{g^) (48) 

Notice that in contrast to the evaluation of the Hamiltonian, explicit expressions for R\ 
and i?2 are required to compute one-loop corrections to operators. In terms of eigenstates 
of Hq these are 

{m\R^\n) = ie{M - |£;^„|)i!^i^ (49) 



and 



{m\R2\n) = /(^i \^-rn\) ((^^\V2\n)+J2{m\V,\q){q\V,\n)QT{E,n,E,m;Ai)] • (50) 



E 

■'-'nm 



17 



Calculating radiative corrections to operators thus becomes a straightforward, albeit tedious, 
process of evaluating commutators. Note that the similarity transformation may be iterated 
to form the renormalization group improved effective operator, just as for the Hamiltonian. 
Matrix elements of an operator O, 

{m\0\n) = {m\S-^OeffS\n) (51) 

have mixings to eigenstates with large energy included perturbatively in Oeff This enables 
one to resolve an ambiguity in applying perturbative QCD to hadronic matrix elements. 
The ambiguity is illustrated in Fig. 1 where two possible diagrams contributing to the 
hadronic matrix element of a generic current are presented. The question is if the current 
couples to nonvalence structure in the hadron (say, the strangeness content of the proton) , 
does it do so via an existing nonvalence component of the wave function (as in Fig. la) 
or via a virtual quark loop (as in Fig. lb)? It is impossible to resolve this in typical 
approaches to this problem {i.e., without reference to dynamics). The method presented 
here allows one to identify both contributions. Fig. la follows from solving the bound state 
effective Hamiltonian in a basis which allows the appropriate Fock state mixing. Thus the 
intermediate states associated with the wave function exist below A. Alternatively, Fig. lb 
arises from transforming the current O — > Ogff and involves intermediate states which exist 
above A. We emphasize that correctly interpreting this situation can only be done in the 
context of a consistent dynamical model. 

FIGURES 

a. 

Fig. 1. Resolving Ambiguities in Current Matrix Elements. 
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A. The Effective Quark Operator 



As a first simple application of the similarity transformation we evaluate the one-loop 
correction to the bare quark operator, 6(k, r). This means evaluating the commutators 
indicated in Eq. with O = 6(k, r). Since we are interested in the evolution of the operator 



itself (and not the generation of new operators such as 6a^), the first order commutator 
([-Ri, 6(k, r)]) does not contribute. Furthermore the order commutator, [_R2, &(k, r)], is 
zero. This is because the first term in R2 (see Eq. [Q) involves the Coulomb potential 
which is normal ordered. Thus all commutators with the quark annihilation operator simply 
yield terms which are proportional to the Coulomb interaction (and never 6(k, r)). The 
commutators of the second term in R2 cancel in pairs and again, do not contribute. Thus 
the only contribution to fog// is the double commutator over Ri in Eq. Two of the 

64 possible terms are nonzero, these are shown diagrammatically in Fig. 2a,b. As Fig. 2 
indicates, the process of evaluating commutators for Oe// is similar to performing time- 
ordered perturbation theory. However, we warn the reader that this is an analogy only since 
nontrivial cut-off dependence is implicit in each diagram. 




i 



a. b. c. 

Fig. 2. Diagrams Contributing to beff- 

Expressions for these contributions to fog// are respectively 

_ ,Cf r dp (l-p-lk-i)e(Ao>|E(p)+^(I)-i^(k)|>Ai) 

^ 4 i (27r)3 ^(1) (E(p)+^(I)-E(k))2 ^ 

^Cf f dp (l + p-ik-i))e(Ao > |S(p)+o;(l) + E(k)| > Ai] 



^ 4 i (27r)3 ^(I) {E{p)+uj{l) + E{k)y ' ^^^^ 

here we have taken 1 = p — k and Q{a > b > c) = 6{b — c)6{a ~ b). Applying the similarity 
transformation ~ log Ao/A^r times to reduce the cut-off from Aq to A^v yields the following 
leading order result 
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,2 



Kk,x)^Kk,x)(^) 



Cfj£ 



(47r)^ 



(54) 



When the similarity scale equals the ultraviolet cutoff, no renormalization of the quark 
operator occurs as expected. As the similarity scale is reduced, the "strength" of the quark 
operator diminishes as "non-diagonal" operators (such as ba'^) become important. 



As a simple application of the renormalization of composite operators, we now apply 
the methodology to the axial charge operator for massless quarks. In terms of bare quark 
operators this is given as 



In this case the first six diagrams of Fig. 3 contribute to the perturbative corrections to 
Q5 (where now the dot refers to O — Qs). The Coulomb loop diagrams (Figs. 3.g and h) are 
zero because does not contain terms proportional to b^d^ or bd. Doing the requisite spin 
sums and performing the integral over loop momentum and k leads to the same integrals 
for Figs. 3. a and 3.b as in Eqs. [52] and [53]. These are cancelled by the vertex correction 
integral of Fig. 3.c. The same holds for the Z-graphs of Figs. 3.d,e, and f. Thus there is 
no correction to at order g^. This is as expected because the axial current is partially 
conserved and carries no anomalous dimension. 



B. The Axial Charge at One-loop 




(55) 



20 




h. 



Fig. 3. Diagrams Contributing to Oeff- 



C. Renormalization of 



The evaluation of Eq. |^ in the case of the qq operator is more involved than in the 
case of the axial charge because diagrams which correspond to b^b and Sd also contribute. 
The regulated operator qq^° is 

dq 



E 



-E 



(27r)= 



s(q)[6t(q, r)6(q, r) + d\-q, r)rf(-q, r)] 



-c{q)e{^ - i?(q))[6t(q, r)dt(-q, r) + d\-q, r)6(q, r)] 



dc[ 



(27r)3 

where s(q) and c(q) are given by 



c(q) 



mo 



g2 + m-Q 



|q| 



g2 + mQ 



(56) 

(57) 
(58) 
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After similarity transformation the effective operator qq^ff can be written as 



q(t = + 9 



{471^^ J (27r)3 
- Ae(k; Ao, Ai) [b\k, T)d\-k, r) + d{-k, T)b{k, r)]] 



A,(k; Ao, Ai) [6t(k, T)6(k, r) + d^{-k, T)d(-k, r)] 



(59) 



Recall that qq^° = qqejf by definition. The integrands in this expression may be separated 
into portions representing contributions from [i?2,??] and respectively: 



^ = ^2 + ^11, ^ = ^2 + ^11- 



(60) 



The commutator qq^^] leads to an effective operator that changes the number of gluons 
by ±1 and hence need not be considered. The parts of these commutators which give rise 
to the leading renormalization group trajectory are given by 



A,2(k;Ao,Ai) = s(k)log 



Ao 



^-W{\k\+E{k)-K^)[l- 



K,-E{ky 



+ s 



+ 



(k)/d|q|d(q-k)^^(^-|q|) 

g(Ao > \E{k) - E(q) + a;(k - q)| > AQ 
E(k)-E(q)+a;(k-q) 



^(Ao > I - E{k) - E{q) + c^(k - q)| > Ai) 



-E(k)-E(q)+aj(k-q) 



(61) 



and 



A2(k;Ao,Ai)= c(k)log^ 



+ /d|qM(q-k)^(^-|q|) 
+ (c(k)-k-q) 



8-4^(|k| + E(k) - Ai) 1 



Ai-£;(k)^ 



(c(k) + k • q) 



^(Ao > I - E{k} - £(q) + ^(k - q)| > Ai) 



-E{k)-E{q)+u;{k-q) 
e(Ao > |£;(k) - £;(q) +a;(k- q)| > AO 



E{k)-E{q)+u;{k-q) 



(62) 



The remaining terms are 



An 



and 



^n(k;Ao,Ai) = -2s(k)log:|^-s(k) J ci|q|ci(q • k)[...] 



yleii(k;Ao,Ai) = -2c(k)log^ - J d|q|d(q • k)[...]. 



(63) 



(64) 
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The two integrands in Agu and Acu abbreviated by [...] are identical with the integrands 
in As2 and Ac2 respectively and therefore cancel in Ag and A^. The contributions to As2 
and Ac2 which are proportional to ^^(|k| + E(k) — Ai) vanish when matrix elements of qqeff 



are taken in a basis in which -ffg// is nonzero i.e., for Ai > 2E(k) (see Eq. [|T9|). The 
additional contribution to R2 resulting from off-diagonal self energies {Gc, Gt), discussed 
below Eq. results in terms proportional to 6{2E{'k) — Ai) and also vanishes below the 
cut-off. Finally, after a single application of the similarity transformation, 



qq 



Ao 



1 + 9 
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Ao 
Ai 



qq^ + "finite" 



1 + 9' 



(47r) 



Nt 



qq^K (65) 



The term "finite" refers to contributions which are finite in the limit Aq 00, with fixed 
Ai, and thus do not contribute to the leading renormalization group trajectory. Applying 
the similarity transformation ~ logAg/Ajv times yields. 



qq 



Ao 



A AT/ 



(66) 



Using Eq. [43] we may easily show that the product mqq, is constant along the renormal- 
ization group trajectory i.e., as Aq 



^„Aa 



m{Ao)qq ° = 'm{A]^)qq 



(67) 



IV. NONPERTURBATIVE PROPERTIES 

It is known that constituent quarks are the appropriate degrees of freedom for the de- 
scription of hadrons. Since the DQM is constructed with current quarks, constituent quarks 
must be dynamically generated as the pseudoparticles of the model. We therefore choose 
to proceed by making a pairing ansatz for the vacuum of the standard BCS form. This 
is motivated by the constituent quark model as we shall see shortly. The incorporation 
of a nontrivial pairing vacuum may be easily achieved by performing a Bogoliubov-Valatin 
canonical transformation on the quark fields: 
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5(k, r) = cos(^0(k))6(k, r) - /i(r) sm{^<jy{k))d\K r) 

D(k, r) = cos(^0(k))ci(k, r) + h{r) sin(^0(k))6t(k, r) (68) 

where /i(r) = ±1 is the hehcity. The transformation is parameterized in terms of the 
Bogohubov angle, 0(k). In future we shall use the notation 

5i(k)=sin(^0(k)), 5(k) = sin(0(k)), (69) 

with similar expressions for the cosine. 

Expressing the effective Hamiltonian, H^^^ in terms of the constituent quark and an- 
tiquark operators {B and D respectively) leads to a one-body off-diagonal operator of the 
type 

r flic r n 

^ J (2^^^^' '^^ [^'^^' ^) + ^^^^^^ • 

The Bogoliubov angle is chosen so that this operator vanishes 

F(k» = 0; (71) 

this is the vacuum gap equation of the DQM. 

Although the gap equation allows the determination of a nontrivial vacuum, and hence 
of the dynamical quark mass (discussed below), there is actually another important reason 
for imposing this structure on the model. In particular, the successes of the constituent 
quark model indicate that valence quarks saturate hadronic states to a large degree. The 
BCS ansatz is one of the simplest methods which allows for dynamical chiral symmetry 
breaking and decoupling of qq pairs from the vacuum. Thus the nonperturbative approach 
we have taken is firmly based on the successes of the CQM. We note that the decoupling 
of the vacuum does not persist as states with four or more constituents are included in the 
analysis. In this case it is only possible to approximately diagonalize the Hamiltonian by 
further changing the basis as it is done, for example in the Random Phase Approximation 
(RPA) or by using bound state perturbation theory with either Tamm-Dancoff or RPA 
eigenstates as a basis. 
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The explicit form of the gap equation is given by: 

Sik) + e{^-Eik))Giik;A)Cik) 



{vij\K q) + V?{K q)) 25(k)^? (q) + {vS\k, q) + v!r'\k, q)) C(k)^(q) 



'k2 + m(A)2 + Si(k,A) 
(iq 

+ (V,?)(k, q) + V!r'\k, q)) S{k)S{q) + (V^)(k, q) + V^')(k, q)) 2C(k)S| (q)] . (72) 

Here the functions V^l^ and V^*'* come from the static, two-body, Coulomb and linear (denoted 
by the subscript cl) part of the effective Hamiltonian and the effective interactions resulting 
from transverse gluon exchange (denoted by T) respectively. They are all summarized in 
the Appendix. The only self-energy type contributions, Hi and Gi, come from the linear 
potential. As discussed in Sec. IID, these are necessary to maintain the IR stability of 
physical observables. The importance of the linear self energy term has been noted previously 
in Refs. P JT^JT3|] . It should be noted that the solution to the gap equation depends on the 
scale A, = 0(k; A). 

Since the gap equation incorporates all order g'^ physics correctly it is completely UV 
finite. This is not true for earlier work. For example, the gap equation of Le Yaouanc et 



al. [|14| may by obtained from this one by setting all theta functions equal to unity and by 
ignoring all transverse gluon terms. This introduces an ultraviolet divergence due to the 
Coulomb self-energy term which they avoid by simply neglecting the Coulomb interaction. 



Adler and Davis [|1^ discuss the gap equation with a pure Coulomb interaction. It should 



be noted, however, that in absence of transverse gluons, pure Coulomb exchange leads to 



momentum dependent counter-terms. Finally, the gap equation of Finger and Mandula |T5 
is similar to those above except that they ignore all self energy terms. This eliminates all 
UV divergences but introduces IR divergences. We stress that the only consistent way to 
derive a gap equation which incorporates gluons is with a well-defined Hamiltonian-based 
renormalization scheme, such as is presented here. 

As |k| ^ oo the solution to the gap equation approaches the perturbative result: S{k) — *• 
s(k). It is thus natural to define an effective quark mass with the relationship 

s{k-mdyn{k)) = S{k). (73) 
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The dispersion relation for such a dynamical quark mass is shown in Fig. 4 (the bare quark 
mass was chosen to satisfy m(l GeV)= 5 MeV). The solid hne is E{k;A) for A = 4 GeV. 
The dashed line represents ^Jk"^ + m2(4GeV) (where the quark mass has been perturbatively 
run from 1 GeV to 4 GeV). One sees that the correct high energy behavior is recovered 
and that a constituent quark mass of roughly 100 MeV is obtained at low energy. This 
is approximately one half of the constituent masses used in relativistic quark models. In 
general the constituent quark mass will tend to zero as the scale is reduced and will saturate 
as the scale is increased. A detailed study of the dynamical quark mass is beyond the scope 
of this paper and will be presented elsewhere. 
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k [GeV] 

Fig. 4. The Dynamical Quark Dispersion Relation 

With a nontrivial solution to the gap equation we may proceed to define the nonpertur- 
bative quark condensate. The matrix element of qQeff in the nonperturbative vacuum, \Q) 
(defined hy B\Q) = D\fl) = 0) is given by (per quark fiavor) 

(nmn) = -6 y ^ [^(k)c(k) + {C{k) - l)s{k)] . (74) 
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This expression has been regulated by subtracting the perturbative contribution: 

r dk 

(0|gg|0) = -ey j^A^)- (75) 

The condensate is shown as a function of the renormahzation scale in Fig. 5. The circles 
represent the solution for massless quarks while the diamonds are for m(lGeV) = 5 MeV. 
The perturbative solution is shown as a solid line (this has been normalized to match the 
nonperturbative results). As expected, both nonperturbative curves approach the pertur- 
bative solution as A increases. The behavior of the condensate below 3 GeV is driven by 
the nonperturbative renormahzation group trajectory of the model. How closely this follows 
the actual nonperturbative evolution of QCD is a function of how good our approximations 
are (in particular the form of Kon/ and the vacuum Ansatz). A way to test the consistency 
of the model is to calculate the nonperturbative running quark mass and strong coupling 
constant (and wave function normalization). This may be done, e.g., by fixing A (below 
3 GeV), calculating several observable (such as the pion and rho masses) and fixing the 
parameters of the model to reproduce experiment, 

m^/p(m(A), a5(A); A) = m^/p|expf (76) 

Once the nonperturbative running mass is determined one may compare it to the expected 
running mass (taking advantage of the renormahzation group invariance of mqq and the 
results of Fig. 5.) to assess the accuracy of the model. 

The value of the condensate at low energy is roughly (—100 MeV)^. This does not 
compare well with the commonly quoted value of (—250 MeV)^ [||. The underestimate for 
both the constituent quark mass and the condensate is most likely due to an inadequate 
model for the vacuum. More quark correlations may need to be incorporated or a coupled 
quark-gluon vacuum Ansatz may be required. Unfortunately the use of nonperturbative 
methods implies that the problem is not restricted to the vacuum sector. For example, the 
Thouless theorem relates the pairing vacuum to the RPA pion via the Gell-Mann-Oakes- 
Renner relationship /^m^ = —2mg{qq) 0]. Thus, if the pion mass is fixed, a poor model of 
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the vacuum will be reflected in a low value of the pion decay constant. How severe these 
problems are will be examined in future work. For now, we must remain satisfied that 
nonperturbative statements about the vacuum can be made at all. 




Fig. 5. The Quark Condensate. 
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V. CONCLUSIONS 



This paper presents the development of the Dynamical Quark Model. The model is 
meant to describe hadronic physics and is therefore designed to be amenable to many- 
body techniques from the outset. The use of field theoretic methods forces one to choose 
appropriate degrees of freedom and many-body approximations. The constituent quark 
model has been used as a guide in this process. For example, the choice of the Coulomb 
gauge, the form of the gap equation, and the use of the Tamm-Dancoff approximation are 
all suggested by the phenomenological success of the CQM. 

There are three main stages in the construction of the DQM. The first is the evalua- 
tion of the renormalization group improved QCD Hamiltonian to order g^. This is achieved 
with a Hamiltonian-based renormalization scheme introduced by Glazek and Wilson (the 
similarity scheme). In the second stage, the RGI Hamiltonian is supplemented with a phe- 
nomenological potential. The fact that the starting point of the DQM is QCD provided some 
important constraints. In particular, the Dirac structure of the phenomenological potential 
is determined in the heavy quark limit. Performing the Bogoliubov-Valatin transformation 
comprised the last stage of this program. This allows the construction of a nontrivial vacuum 

- a necessary feature of the model since constituent quarks must be dynamically generated. 
A beneficial corollary is the appearance of chiral pions. 

Since this paper is meant to describe the perturbative development of the model, only 
the simplest preliminary applications have been discussed. Thus we have derived and solved 
the vacuum gap equation, used the resulting Bogoliubov angle to determine the dynamical 
quark mass, and evaluated the chiral quark condensate. It is satisfying and nontrivial 
that a constituent quark mass arises and is roughly as expected from phenomenological 
considerations. This is an indication that the approach we are taking may be useful for 
light quarks (the DQM may be trusted for heavy quarks since it reduces to the CQM in this 
case). Unfortunately, the chiral condensate and the constituent quark mass are too small 

- evidence that the nonperturbative structure of QCD is not fully captured in our model 
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and approximations. Whether this is a serious problem or not can only be determined by a 
comparison with more observables. The application (and possible modification) of the DQM 
to hadronic observables is clearly an important task for the future. 

The DQM has a rich structure and therefore allows the examination of many topical issues 
in hadronic physics. There are several qualitative aspects of the model which we hope to 
explore in the near future. These include an extension of the phenomenological potential to 
incorporate nonperturbative gluonic fiux tube degrees of freedom and an analysis of hadronic 
decays in a spirit similar to our discussion of the dynamical generation of effective scalar 
spin-dependent interactions §]. It will also be interesting to examine hyperfine splittings in 
the N-A and vr-p systems as a function of quark mass. This should shed some light on the 
genesis and utility of the constituent quark model. 

Finally, there is a wealth of specific phenomenological problems in hadronic physics 
which the DQM may address. For example, since gluonic degrees of freedom are explicitly 
included in the model, we may examine the properties of glueballs and hybrids - including 
their couplings to quarkonia. Pions are ubiquitous in hadronic physics, appearing as strongly 
interacting probes, decay products, and exchange currents. It is therefore essential that they 
be thoroughly understood. However the chiral and relativistic nature of pions have made 
this a longstanding deficiency of microscopic models of QCD. It will therefore be interesting 
to examine a selection of problems which involve pions (recall that pions are relativistic 
pseudo-Goldstone bosons built from quark pseudoparticles in the DQM). In particular we 
plan to calculate the pion mass, decay constant, electromagnetic form factor, and the width 
for ttq — > 77. Many other problems may be addressed with the DQM. We regard the scheme 
presented here as the next step in a process where models are used to elucidate and guide 
experiment which in turn directs the development of newer and more sophisticated models 
of low energy QCD. 
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APPENDIX A: 



The functions Sc and Gc appearing in Eq. |115| are given by 



Sc(k; Ao^) = g'Cp I -0^,O{A^ " E{<D) (.(k)s(q) + k ■ qc(k)c(q)) K(k, q) (Al) 



Gc{k- A?) = -/Cp J 7^^(Ao^ - E{cD) (c(k)s(q) - k . qs(k)c(q)) K(k, q) (A2) 



(2vr) 



where 



mq) = ^^^3^- (A3) 

The self energies due to transverse gluon exchange are given by 

ST(k; Ao, Al) = -g'Cp J ^ (l - s{k)s{q) - k ■ Iq ■ ic(k)c(q) 
^(Ao>|E(k)-i^(q)-^(l)|>AO _ _ ^ _ J J 

e(Ao>|i?(k) + E(q)+^(l)|>Ai) 



2a;(l) (cu(l) + E(q) + ^(k)) 



(A4) 



and 



Grik; Ao, Al) = g'Gp J -0^-L.e{Ao - \E{c[) + E{k) + a;(l)|)^^(Ao - \E{ci) - E{k) + c.(l)|) 
X (c(k)5(q) - k • iq ■ ic(q)s(k)) eT(^(q) + ^(k) + E{q) - E{k) + u;(l); Ai). 

(A5) 
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The self energies due to the hnear potential appearing in Eq. and Gi, are obtained 

from Sc and respectively by substitution 

47rcr 



^7^Ci.V;(k,q)^l^(k,q) 

|K q| 

where cr = 0.18 GeV^ is the string tension. 

The functions V^;^ which appear in Eq. [^] are given by 

Vi/^(k,q) 
Vi?(k,q) 

Vi?(k,q) 
Vif^(k,q) 



(s(k)s(q)+k-qc(k)s(q))KKk,q) 

^ (c(k)c(q) + k ■ q(l + 5(k).(q))) e_(k, q; A)K/(k, q) 
+ \ (c(k)c(q) - k ■ q(l - .(k)s(q))) 0+(k, q; A)KKk, q) 



-s(k)c(q) + k . qs(q)c(k)J 0q(q; A)KKk, q 
Vi?(q,k) 



(A6) 



(A7) 



and 



V^'^(k,q) 



4cu(l) 



1 - s(k)s(q) - k • iq ■ ic(k)c(q) 
X [ 0T(i?(q) - i?(k) + ^(1), E(q) - E(k) + ^(1); A) 

+ eT(-^(q) + ^(k) + cu(l), -E(q) + E(k) + uj{\)- A)] 
[l + s(k)s(q) + k • iq • ic(k)c(q) 
X [ 0T(-i5(q) - i5(k) + 07(1), -E(q) - E(k) + ^(1); A) 
+ er(S(q) + ^(k) + cu(l), E(q) + ^(k) + uj{\)- A)] 



(AS) 



V^'H^q) 



2uj{\) 



c(k)c(q) - k ■ Iq ■ 1(1 - s(k).(q))J e_(k, q; A) 

X eT(-^(q) + ^(k) + cu(l), E(q) - ^(k) + cu(l); A) 

c(k)c(q) + k ■ iq ■ i(l + s(k)s(q))] 0+(k, q; A) 
X eT(^(q) + ^(k) + u;(l), -E(q) - E(k) + cu(l); A) 
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Vi?)(k,q) 



.(k)c(q)-k-iq-ic(k)s(q))] e^{q,A) 
X [ Ori-Eici) - Eik) + coil), i5;(q) - Eik) + co{iy, A) 
+ eT(^(q) + ^(k) + ivil), -E{q) + E{k) + A)] 



V^f)(k,q)=V^^^(q,k) 



(3), 



where Vd is the sum of Coulomb and hnear potentials 

KKk,q) = + 



(A9) 



(AlO) 



2|k-q|2 |k-q|4 

and the effective transverse gluon cut-off Gt is given by Eq. [^]. Additional cut-offs are 
given by 



e±(k,q;A) = 0(A-2|£;(k)±E(q)|) 



(All) 



and 



Qp = 9{A-2E{p)). 



(A12) 



Furthermore s(k) and c(k) stand for the sine and cosine of the free particle BV angle i.e. 

1 



s(p) = sin(0„(p)) 



m 



C(p) = COs(0m(p)) 



(A13) 



Finally, 1 = k — q and a;(l) = |1| represent the momentum and energy of an exchanged gluon 
respectively. 
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Fig. 1. Resolving Ambiguities in Current Matrix Elements. 
Fig. 2. Diagrams Contributing to ftg//- 
Fig .3. Diagrams Contributing to Oe//- 

Fig. 4. The Dynamical Quark Dispersion Relation. The sohd hne is the numerical result 
while the dashed line is the perturbative relation for m{ 4 GeV) = 3.1 MeV. 

Fig. 5. The Chiral Condensate. The sohd hne is perturbation theory, the symbols 
represent the nonperturbative condensate for massless quarks (circles) and for m(l GeV) = 
5 MeV (diamonds). 
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